Wave scattering by discrete breathers 
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We present a theoretical study of linear wave scattering in one- dimensional nonlinear lattices by 
intrinsic spatially localized dynamic excitations or discrete breathers. These states appear in various 
nonlinear systems and present a time-periodic localized scattering potential for plane waves. We 
consider the case of elastic one-channel scattering, when the frequencies of incoming and transmitted 
waves coincide, but the breather provides with additional spatially localized ac channels whose 
presence may lead to various interference patterns. The dependence of the transmission coefficient 
on the wave number q and the breather frequency Qt is studied for different types of breathers: 
acoustic and optical breathers, and rotobreathers. We identify several typical scattering setups where 
the internal time dependence of the breather is of crucial importance for the observed transmission 
properties. 
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I. INTRODUCTION 

The problem of wave propagation through media with 
various inhomogeneities has been a complex issue of con- 
stant interest and appears in different areas of physics. 
Particular examples are acoustic and electromagnetic 
wave propagation in various disordered media [1,2], tun- 
neling of electrons in solids [3] and electron transport 
through quantum (molecular) wires [4,5]. In many cases 
of interest the conductivity (electron transport) [3,5,6] 
and the heat conductivity (phonon transport) [7-9] are 
determined through the wave scattering by spatial inho- 
mogeneities. Of particular interest is the wave propaga- 
tion in one-dimensional systems where interference effects 
may be strongly enhanced. 

In most of the studies wave scattering by static local- 
ized inhomogeneities has been considered. More recently 
the scattering by generic time-dependent potentials has 
received strong attention [10-12]. This is due to the pos- 
sibility to generate various time-dependent scattering po- 
tentials artificially e.g. in the presence of laser beams or 
microwave radiation. Several interesting effects such as 
giant enhancement of tunneling [10,11] and Fano reso- 
nances [13-16] have been found. 

It is a well established fact that various nonlinear spa- 
tially discrete systems can support different types of exci- 
tations, namely, propagating linear waves (phonons) and 
time-periodic spatially localized excitations called dis- 
crete breather states (DB) [17-19]. The origin of the lat- 
ter localized states is not the presence of disorder but 
rather the peculiar interplay between the nonlinearity 
and discreteness of the lattice. While the nonlinearity 



provides with an amplitude-dependent tunability of oscil- 
lation or rotation frequencies of DBs, Qt, the spatial dis- 
creteness of the system leads to finite upper bounds of the 
frequency spectrum of small amplitude waves LUq. It al- 
lows to escape resonances of all multiples of the breather 
frequency Q(, with LUq. Note here, that nonlinear dis- 
crete lattices admit different types of DBs depending on 
the spectrum of linear waves propagating in the lattice, 
i. e. acoustic breathers and rotobreathers (Fig. la and 
lb), optical breathers (Fig. Ic), etc. Such properties of 
DBs as their frequency dependent localization length and 
the stability of DBs with respect to small amplitude per- 
turbations have been widely studied. DBs have been ob- 
served in experiments covering such diverse fields as non- 
linear optics [20], interacting Josephson junctions [21,22] 
magnetic systems [23] and lattice dynamics of crystals 
[24]. 
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FIG. 1. A schematic representation of different types of 
discrete breathers: a) acoustic breather; b) acoustic roto- 
breather; c) optical breather. 

Although DBs present complex dynamical objects, in 
many cases experimental measurements can be well un- 
derstood by using certain time-averaged properties of 
DBs [25,26]. Thus, a natural question appears whether 
the internal breather dynamics is of crucial importance to 
understand the lattice dynamics in the presence of DBs. 
In this paper we study the propagation of small ampli- 
tude plane waves in one-dimensional nonlinear lattices in 
the presence of a DB and obtain that the internal dynam- 
ics of the DB may lead to a drastical increase or decrease 
of wave transmission as compared to the time-averaged 
setup. Thus the wave scattering by DBs is interesting 
both as some spectroscopical tool to study DB properties 
and as a way to control the wave transmission (conduc- 
tivity) by varying the DB state. Finally our studies are 
of use for the general understanding of wave scattering 
by time-periodic potentials. 

First successful attempts to describe the variety of phe- 
nomena arising from wave scattering by DBs have been 
performed some time ago [27 30,16]. While a number 
of results have been obtained, we are far from a full 
description of the complexity of possible scattering out- 
comes. This concerns the wave scattering in systems with 
acoustic spectra uig, the importance of the internal DB 
dynamics including the comparison between rotational 
and vibrational excitations, the dependence of the wave 
propagation on the DB energy, and the single channel 
elastic versus the two channel inelastic scattering cases 
[28]. Here we will address these problems but restrict 



ourselves to the elastic scattering case. 

The paper is organized as follows: in Section II we 
present the general formalism and describe analytical and 
novel numerical methods used to analyze the wave scat- 
tering by DBs, in Sections III, IV and V the dependencies 
of the transmission coefficient on the breather frequency 
fif, and the wave number q for different types of breathers 
(acoustic breathers and rotobreathers, optical breathers, 
see Fig. 1) are obtained, and finally a discussion is pro- 
vided in Section VI. 

II. GENERAL FORMALISM 

To proceed we will consider one-dimensional nonlinear 
lattices with nearest neighbor interaction, optional on- 
site (substrate) potential, and with one degree of freedom 
per lattice site. Both the increasing of the interaction 
range and the extension to more than the one degree of 
freedom per lattice site are not of crucial importance. 

The dynamics of the system is characterizcid by time- 
dependent coordinates Xn{t) and the class of Hamiltoni- 
ans considered here reads 

^^J2f'A + V[Xn] + W[Xn-Xr,-l]] . (1) 
n ^ ' 

Here V"[.t] is an optional on-site (substrate) potential and 
M^[x] is the nearest neighbor interaction. The equations 
of motion become 

Xn = -W'[Xn - Xn-l] + W'[Xn+l - Xn] - V'[Xn] ■ (2) 

Without loss of generality we take V[0] = W[0] = V'[0\ = 

W'[0] = and V"[0] > 0, W"[0] > 0. This Hamiltonian 
supports the excitation of small amplitude linear waves 
with the frequency spectrum 

u^l = V"[0]+m"[0]sm'{£) , (3) 

with q being the wave number. 

Time-dependent spatially localized solutions (DBs) ex- 
ist for different types of potentials V[x] and al- 
though at least one of the two functions {V[x] and/or 
W^[x]) has to be anharmonic. DB solutions are charac- 
terized by being time-periodic Xn{t+Ti,) = Xn{t) and spa- 
tially localized a;|„|^oc ^ (except systems with V = 
where Xn^±oc ±f' with c possibly being nonzero). If 
the Hamiltonian H is invariant under a finite translation 
(rotation) of any Xn ^ a:„-|-A, then discrete rotobreathers 
(DRB) may exist [31]. These excitations are character- 
ized by one or several sites in the breather center evolv- 
ing in a rotational state xo{t + T^) = xo{t) + A, while 
outside this center the lattice is governed again by time 
periodic spatially localized oscillations. The breather fre- 
quency ilfc = 27r/Tf, can generally take any values pro- 
vided k^lb 7^ Uq for all nonzero integer k. As is an 
analytic function of q, DBs are exponentially localized 
on the lattice. 



2 



A. The linecirized problem 

To study the scattering of small amplitude plane waves 
by a DB wc linearize the equations of motion (2) around 
a breather solution a;„(t) = x„(t) + e„(t): 

'in = -W"[xnit) - x„_i(t)](e„ - e„_i) 

+W"[Xn+l{t) - £nmen+l - £„) - V"[Xn{t)]en ■ (4) 

This is a set of coupled linear differential equations with 
time periodic coefficients of period T^. Note that these 
coefficients are determined by the given DB solution 

Xn{t). 

Eq.(4) determines the linear stability of the breather 
through the spectral properties of the Floquet matrix, 
[19,28] which is given by a map over one breather period 

where e = (..., e„, e„+i, ...). For marginally sta- 
ble breathers all eigenvalues of the symplectic matrix F 
will be located on the unit circle and can be written as 
e*^. The corresponding eigenvectors are the solutions of 
Eq.(4), and fulfill the Bloch condition 

en{t + n)=e-'''/^^en{t). (6) 

Because the DB solution is exponentially localized on the 
lattice, equation (4) will reduce to the usual small am- 
plitude wave equation far away from the breather cen- 
ter. Thus, only a finite number of Floquet eigenvectors 
are spatially localized, while an infinite number of them 
(for an infinite lattice) are delocalized and these solu- 
tions correspond to plane waves with the spectrum (3) 
and the Floquet phases 9 = ±cOqTi,. The remaining Flo- 
quet eigenvalues correspond to local Floquet eigenvectors 
and are separated from the plane wave spectrum on the 
unit circle. Note here, that two eigenvectors with the de- 
generated eigenvalue e'^ = -|- 1 always exist and reflect 
perturbations tangent to the breather family of solutions. 

As a consequence of the Bloch condition (6) any spa- 
tially extended Floquet eigenvector can be represented in 
the form 

oo 

en{t)= Yl e„fee'(-'+'="'')*. (7) 

fc=— oo 

What happens if we send a plane wave with frequency 
u!q to the DB? We will deal with the case of one-channel 
scattering as for any k ^ and any q' , uiq' ^ Uq + kilb- 
This condition determines that all channels with nonzero 
k are 'closed', i.e. the spatial amplitudes tnk are localized 
in space. Note here, that the frequencies of transmitted 
and reflected waves are the same and coincide with the 
the frequency of the incoming wave, since they all belong 
to the only open channel with A; = 0. 



For harmonic interaction potentials W it was shown in 
Ref. [28] that the momentum 

J = -PF"[0](Ime;e„_i) (8) 

is independent of the lattice site. Here the averaging 
is meant with respect to time. In a similar way it is 
straightforward to show that the quantity 

J = -{W"[xn{t) - x„_i(t)]Ime;e„_i) (9) 

is independent of the lattice site. Since the breather so- 
lution Xn(t) is spatially localized, at large distances from 
the breather we find again that the momentum J is in- 
dependent of the lattice site. Especially we find that it is 
the same for n — > ±oo. Following Ref. [28] we conclude 
that the one-channel scattering case is elastic, i.e. the en- 
ergy of an incoming wave equals the sum of the energies 
of the outgoing (refiected and transmitted) waves. 

Despite the fact that we will study the one-channel 
elastic scattering, the scattering potential of DB is time- 
periodic. The above mentioned 'closed' channels are ac- 
tive inside the breather core, i. e. the frequency of linear 
waves in a finite area around the breather center can 
change due to the interaction with the DB (see Fig. 2). 
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FIG. 2. Schematic representation of the one-channel scat- 
tering of a wave by a discrete breather. 

Thus one of the questions to be answered below is to 

identify the cases when the well-known scattering by a 
time- averaged (static) potential is not sufficient to de- 
scribe the actual results of wave scattering by discrete 
breathers. This implies that interference effects through 
local interactions between the active closed channels may 
substantially change the scattering results as compared 
to a time-averaged scattering potential. 

Assuming that the breather is located around the cen- 
tral site n = 0, the one-channel scattering problem can 
be written as 
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e„(t) = ^re-*^'"'*-^") ,n > (10) 



for 



|n| >sup[6(0),6(a;,)] 



(11) 



where Cb(u;) = sup^(a; + kilb) and 



sinh2eM/2 = ^g2Lj±! , |H<o;,(0) 
cosh^^M/2 = 4^ , W\>iv,{n) 

^ measures the characteristic inverse locahzation length 
of a closed channel at frequency v (note that for u = kflj, 
the localization length is that of the breather itself, see 
also [28]). The incoming wave has amplitude A[, and 
the reflected and transmitted wave amplitudes are given 
by Af{ and At respectively. The transmission coefficient 

tg = \AT/Al\\ 

For further considerations we will use the notation of 
the Bloch functions Cn{t) defined as 



<^nit) = C,n{t)e~ 



(13) 



In order to estimate the relative strength of closed 
channels with fc ^ we expand the time-periodic co- 
efficients of (4) in a Fourier series with respect to time: 

oo 

W"[Xn{t)-Xn-l{t)]= Wn,ke''''^\ (14) 

k= — oo 

oo 

V"[Xn{t)]= ^n,fce"="^*. (15) 



fe=— oo 

We consider first the case of a strongly localized optical 

breather located at site n = Q with W{y) — Taking 
into account a single closed channel with some value of 
k, inserting (14), (15) and (7) into (4) and excluding e„fc, 
the relative strength Sk of the closed channel to the open 
one will be given by 



"o.k 



{vo,o - l)(t^o,o + 2c(l - r]k) - {ojq + kflb)"^) 



where the relative amplitude 

Vk 



eo,fe 



(16) 



(17) 



has positive sign for ujq + kCli, located inside the phonon 
gap and negative sign otherwise. For <C 1 we do not 
expect any significant contribution from the given closed 
channel, while Sk > 1 indicate a strong influence of the 
closed channel on the scattering process. Note that the 
expression ^yvoo'+'^c{^^^rikj in the denominator of (16) 



is just a frequency ujl of a local phonon mode of the time- 
averaged scattering potential. For spatially discrete sys- 
tems these local phonon modes may be located inside 
or outside of the phonon gap. The denominator of (16) 
may vanish for certain wave numbers, which would im- 
ply a resonance-like enhancement of the closed channel 
contribution (for certain wavcnumbcrs q) . As such a res- 
onant enhancement of a closed channel amplitude acts 
as a huge effective scattering potential to the open chan- 
nel, for these cases we expect a resonant suppression of 
transmission. Thus, qualitatively such a complete sup- 
pression of transmission (Fano-like resonance [13,15]) is 
explained by a resonant interaction of the propagating 
phonon with the specific local phonon mode. However 
in order to quantitatively analyze this effect the renor- 
malization of the value of ujl due to all nonresonant pro- 
cesses, has to be taken into account. It will be done below 
for a particular case of optical breathers by making use 
of the Green function method. We obtain that although 
the renormalization oftOL is rather small it may become 
important especially as the width of a phonon band is 
small, W" < 1. 

In a similar way we proceed for the estimation of the 
closed channel contribution of acoustic breathers. We 
obtain the following expression for the relative strength 

Tk 



Tk 



%,k 



{wo,o){wo,o - (Wg -f- fcf2b)2) 



(18) 



Thus, we again find a resonant suppression of transmis- 
sion when the presence of acoustic DB leads to a local 

increase in the nearest neighbor interaction and to an 
appearance of a corresponding local phonon mode. How- 
ever, this case is much more involved as compared to 
the optical DB case, and Eq. (18) may serve only as 
a qualitative tool to check whether a closed channel is 
strongly contributing to the transmission or not. We 
will instead provide with a more quantitative analysis 
for the cases considered, based on the particular acoustic 
breather properties. 



B. Numerical scheme 

To compute numerically the transmission coefficient we 
have developed a scheme which generalizes the one given 
in Ref. [28] (which relies on a spatial reflection symmetry 
of the breather and thus of the scattering potential). At 
variance with Ref. [28] our scheme is capable of dealing 
with any (perhaps spatially nonsymmetric) time-periodic 
scattering potential. 

We look for solutions of Eq.(4) which correspond to 
zeroes of the operator 



G(e-I0),i'(0)) = 



m 
m 



(19) 
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on a lattice with 27V + 1 sites labeled -N,{-N + 
1), — 1, 0, 1, ...(iV — 1),N. The incoming wave is fed 
from the left, and the transmitted wave is leaving the 
system to the right. The boundary condition at the right 
end is eN+i = e"**^'*, which implies that the transmit- 
ted wave will have amplitude 1. With a given boundary 
condition at the left end e-at-i = (A + ii?)e~*"«*, where 
A and B are real numbers, we may find the zeroes of 
(19) using a standard Newton routine. Due to the lin- 
earity of the equations of motion in e an arbitrary initial 
guess and one Newton step are needed to converge to the 
zeroes. In practice due to roundoff errors an additional 
Newton step may be required. 

However with arbitrary A and B we will not realize 
the scattering case (10,10) in general. This is due to the 
fact that all extended Floquet states of an infinite sys- 
tem are twofold degenerated because time reversal sym- 
metry holds far from the breather center. To succeed we 
add a second Newton procedure which uses A and B as 
free parameters such that the solution on site N becomes 
ejv — e"*'""^"*, ensuring that we realize a single trans- 
mitted traveling wave of amplitude one at the right end of 
the system. After the Newton procedures are completed, 
the transmission coefficient is then given by 



to 



4 siu q 



\{A + iB)e-^i - C- 



(20) 



While the Bloch functions arc in general time- 
dependent close to the breather center, they will be time- 
independent complex numbers at large distance from the 
breather. We note that the computation operates at the 
machine precision, and we obtain results which are size 
independent, i.e. with the above described boundary 
conditions we emulate an infinite system. The discrete 
breather solution itself has to be obtained beforehand 
using standard numerical procedures [32]. 

The enumerator in (20) vanishes at the extremal values 
of u)q, i.e. at g = and q = n. If the denominator is finite 
at these values of q, the transmission will also vanish. 
This is indeed the generic case for a quadratic dependence 
of the spectrum ojg on q around these points. Exceptions 
are expected for acoustic chains (see next Section). 

Another peculiar point is that if upon changing some 
control parameter, e.g. the breather frequency, a local- 
ized Floquet eigenstate attaches to or disattaches from 
the extended Floquet spectrum, the transmission coefii- 
cient will be exactly i = 1 for the g- value which corre- 
sponds to the edge of the spectrum Wg, i.e. g = or 
q = Tr [28,29]. 



This method is especially convenient as the scattering 
potential is localized in space. To apply this method to 
a particular case when the presence of DBs leads to the 
appearance of a localized time-dependent on-site scatter- 
ing potential, F"[i;„(t)] — 1, we perform a time Fourier 
transformation of Eq. (4) and obtain the equation for 
the Green function Gi^g{ni,n2) ■ 

-Y.GlM,m) j dWn{m)G^^+n{m,n2) , (21) 

where G^,^ (77.1,77.2) is the Green function of the linear 
equations of motion in the absence of the DB. The 
Fourier transform of the localized scattering potential 
Un{m) = J dte'^\V"[xm{t)] - 1) is determined by the 
properties of the DB solution. Because the DBs are peri- 
odic solutions in time the potential Uq (m) contains just 
the harmonics of the breather frequency, fl = kflf,. 
Moreover, as we will see later in the considered cases 
we can take into account the harmonics with small val- 
ues of A; =0 ,±1, ± 2 only. Thus, the calculation 
of the Green function G^^ (ni, 722) can be represented in 
a diagrammatic form where terms which correspond to 
active closed channels describe the local (virtual) absorp- 
tion and emission of phonons by the propagating phonon 
in the presence of the DB (some of the typical diagrams 
are shown in Fig. 3). Moreover, for one-channel scatter- 
ing the energy conservation law of absorbed and emitted 
phonons has to hold. In order to obtain the transmission 
coefficient tg we need also an expression for the Green 
function (721,112) which reads 



G° (721, 712) = - J 



dq e 



ig(ni-Ti2) 



27r a;2-a;2 



(22) 



The Green function of the full problem (for one-channel 
scattering) has a similar form for large distances from the 
breather center (m — > —00 and 722 -f-00): 



Go,, (721,722) = -iDg 



gig\ni-n2\ 

d{u^l)/dq ' 



(23) 



where tq 

As an example for a static scattering potential which 
is strongly localized on a single site {n = 0) we obtain 



Gw, (721,712) = 



l+/3oG" (0,0) ' 



(24) 



C. Green function method for a time-periodic 
localized scattering potential 

We will also analyze the wave propagation through 
DBs by making use of the Green function technique [33] . 



where /3o is the strength of the static potential. The 
transmission coefiicient ig in this case is given by 



{dH)/dqr 
Pl + {d{ujl)ldqY ■ 



(25) 
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We will use the Green function method in Section V 
where the wave scattering by optical DBs (see Fig. Ic) 
will be presented. We show how most important dia- 
grams can be taken into account in the case of a time- 
periodic localized scattering potential. 



with 



Cn- 



( , gn+Cn,n-l-C^, Cn,n-1 



1 







(27) 



(28) 



"1 <H, "0 

I I 



2S^ 



b) 



c) 



no 



fig no "-2 




d) 



FIG. 3. Typical diagrams describing the interaction of the 
propagating phonons with the time-periodic DB scattering 
potential: a) scattering by a time-average (static) potential; 
b) resonant scattering process; c) a process loading to the 
renormalization of a phonon local mode frequency; d) various 
complex processes. Here, thin and thick solid lines present 
the Green functions in the absence of the DB and in the pres- 
ence of the time-average part of the DB scattering potential. 
The tirnc-avorago part of the DB scattering potential is shown 
by a black circle, and the dashed lines show the absorption 
(emission) of phonons. 



D. Time-averaged scattering potential 

Since we are interested in understanding the impor- 
tance of the time dependence of the scattering potential, 
we will also compare the numerical results with those 
obtained by time-averaging the DB scattering potential. 
This time-averaging can be found numerically beforehand 
and Eq. (4) becomes 



in = -WnfiiCn - En-l) + 

Wn+lfii^n+l — En) — Vnfi€n- 



(26) 



Because in this case all inhomogeneities are time- 
independent we can use the standard scattering matrix 
method. With e„(f) = CnC"*'^"* Eq.(26) is rewritten as 



where En = Wn,o and c„,„_i = w„,o- It follows 

C-N 



( Cn+1 
VCjv 



M 



C- 



JV-1 



with 



-N 



M = JJ Mi. 



(29) 



(30) 



i=N 



The expression for the transmission coefficient tq for 
the time averaged scattering potential is then given by 



to 



4 sin q 



\Mri{q)e-ii + Mi2{q) - M2i{q) - M22(g)e'9|2 



(31) 



where M^- are the four matrix elements of the 2x2 matrix 
M. 



III. SCATTERING BY ACOUSTIC BREATHERS 

In this Section we study the wave scattering by so- 
called acoustic breathers. The corresponding systems are 
characterized by a gapless spectrum cjq of propagating 
linear waves, and by a conservation of total mechanical 
momentum. 

The generic choice for the potentials in the Eq. (2) is 



V{y) = and Wiy) 



y + jU ■ This choice 



leads to the well known Fcrmi-Pasta-Ulam system [34]. 
We will first consider the case /3 = which implies the 
presence of a parity in the interaction potential W and 
commcint on the influence of parity violation for (3^0 
later on. The dispersion relation of phonons is given by 



2 sin ■ 



(32) 



A breather solution with frequency Qf, = 4.5 is shown 
in Fig. 4. 
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FIG. 4. Displacements of an acoustic breather with zero 

velocities at t = and Q,h = 4.5. 

Inset: Relative strength ru for the second and fourth closed 
channels versus q. 

The corresponding Fourier components of the scatter- 
ing potential are plotted in Fig. 5. 




FIG. 5. Fourier components Wnk for different k versus n 
for the breather in Fig.4. 

Acoustic breathers can be found with frequencies 
above the phonon band |r^b| > max \ujq\. However, here 
we consider \^lb\ > 2max|Li;g| which implies 1 > 4. 
This condition is needed to realize the one-channel scat- 
tering case of linear waves. For (3 = the ac part of the 
interaction potential W has even harmonics only, and the 
frequency of a propagating phonon in a closed channel 
2f2f, ± uiq can not match the frequency of a local mode of 
the timc-avcraged scattering potential. Thus the closed 
channels play no important role in the scattering pro- 
cess. As a consequence the wave scattering by acoustic 
breathers is practically identical with the scattering by the 
time- averaged potential. Indeed numerical computations 
do not show any relevant difference between the linear 



wave propagation in the presence of time-dependent and 
the static (time-averaged) potentials. Yet there is a num- 
ber of interesting features of the scattering which deserve 
to be exploited. 

For V = the lattice conserves the total mechanical 
momentum P = J2n in addition to the energy. With- 
out loss of generality wc will choose P = here. From 
this conservation law it follows that the total mechanical 
momentum of the linearized problem 11 = J2n is also 
conserved. This implies that for any solution e„(t) the 
shift e{t) -f C is also a solution. In particular e„ = C is 
a solution, which corresponds to a wave with q = and 
for which we have tq^Q = 1. Thus we conclude that the 
transmission of waves by a breather for acoustic systems 
at g = is always perfect, no reflection occurs. This is in 
sharp contrast to systems without conservation of total 
mechanical momentum (see below). 

The peculiar dependence of the transmission coeflicient 
tq on the wave number q and the breather frequency Qi, 
is shown in Fig. 6. The breather frequency is varying 
over nearly three decades. We indeed observe perfect 
transmission at g = 0, and zero transmission at q = tt. 
However we also find that in the studied breather fre- 
quency range two rather narrow peaks around q = n 
appear with perfect transmission. These structures are 
due to the detachment of localized Floquet eigenvectors 
from the continuum of extended Floquet eigenstates. A 
surprising result is shown in Fig. 7. Wc plot the trans- 
mission tq at g = 7r/4 as a function of the breather fre- 
quency We obtain plateaus and crossovers at certain 
breather frequencies. The crossover positions clearly cor- 
relate with the appearance of localized Floquet states, 
which are traced through the perfect transmission close 
to g = TT in Fig. 6. The dependence of the transmission 
on q for the two observed plateaus is shown in the inset 
of Fig. 7. The plateaus range over several decades in Qf,, 
and the g-dependence of tg is rather similar on different 
plateaus. 

To understand this feature wc remind that large 
breather frequencies imply that the breather energy and 
its amplitude in the breather center are large as well. 
Thus we may neglect the harmonic part of the interac- 
tion potential W inside the breather core. The result- 
ing amplitude distribution of the breather core is char- 
acterized by a superexponcntial decay Xn{t) = AnG{t), 
An ^ A^_i, where G{t) is an oscillatory master function 
[35,36]. Due to the symmetry of the breather solution 
the spatial profile is described by j4_„ = —An+i , n < 0. 
The amplitudes An have been computed: Ai — 1, 
A2 ~ -0.166,^3 - 4.796 10--*, yl4 ~ -1.15 10"", 
and so on [34] . The overall amplitude of the oscillations 
Xn{t) is tuned by the amplitude of the master function 
G{t). 

The scattering is essentially described by the time- 
averaged scattering potential. This potential corresponds 
to a large increase in the nearest neighbour coupling 
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terms. In the breather center we find 



0.6 



{W"[xi{t) - ctom = 1 + m^iit) - xo{t)f) (33) 
with 

3{{x,{t) - xoitjf) = ^ (34) 

(the details of the calculation are presented in Appendix 
A). Due to the large value of X for large it becomes 
evident that time-dependent corrections to the scattering 
potential are negligible. By making use of the distribu- 
tion of oscillation amplitudes in the lattice we then find 

3{{Xr^{t) - Xn-l{t)f) ^ AlX\ (35) 

The scattering on a plateau in Fig. 7 is due to a finite 

number of matrices Mn which should be taken into ac- 
count in Eq. (30). Qualitatively the crossover threshold 
can be defined as 

AlX'c^l . (36) 

Thus we obtain the two crossover frequencies flbi = 13.1 
and flb2 = 4527, which match with the observed crossover 
positions. For A < Qi, < ^bi we need only 4 matrices 
Mn, for Qbi < < ^b2, 6 matrices ]VI„ and so on. 
We also computed the transmission at q = 7r/4 using 
the corresponding reduced set of relevant matrices. The 
obtained limiting values of tq are shown as a dashed line 
in Fig. 7. We observe very good agreement between the 
predicted plateau heights and the actual data obtained 
in direct numerical simulations. 




FIG. 6. The dependence of the transmission coefHcient on 
q and Ot. 




1 2 3 4 

10 10 _ 10 10 

ilb 

FIG. 7. The dependence of the transmission coefficient on 

Clt for a particular value of g = ^. The dashed lines show 
the predicted plateau heights (see text). The inset shows the 
dependence of the transmission on q for the two observed 
plateau regions. 

With the above said it becomes also transparent, why 
we observe detachment of localized Floquet states in the 
crossover region. The acoustic breather presents an ef- 
fective potential well for the propagating phonons. The 
width of this well increases with the breather frequency, 
and the number of possible localized Floquet states also 
increases. Such a periodic appearance of perfect trans- 
mission is similar to the quantum mechanical scattering 
by a potential well in the presence of quasi-discrete levels 
[37]. 

Finally we checked the influence oi (3 ^ 0. With re- 
spect to the breather the presence of cubic terms in the 
interaction potential leads to the generation of a kink- 
shaped dc lattice distortion. The corresponding scat- 
tering potential becomes asymmetric in space and odd 
harmonics in the ac scattering potential appear. This 
immediately leads to the possibility of matching between 
the frequency of a propagating phonon in a closed chan- 
nel fib — LOq and several local modes of the time-averaged 
scattering potential. Thus, a resonant suppression of a 
transmission appears in this case. Indeed, we observe 
this effect by direct numerical simulations as shown in 
Fig. 8, where two transmission zeros are found. We will 
discuss this effect in more detail for the case of optical 
breathers below. 
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FIG. 8. The dependence of the transmission coefficient on q 
for an acoustic breather with fit — 4.5 and f3 — 1. The dashed 
line is the result for the time-averaged scattering potential. 
Two transmission zeros are observed at corresponding values 
of q. 

IV. SCATTERING BY ACOUSTIC 
ROTOBREATHERS 

Remarkable differences to the results of the preceding 
Section are obtained if the interaction potential W is cho- 
sen to be a periodic one (the potential V is still zero in 
this Section): 

W{y) = 1 - cos(y). (37) 

With such an interaction potential the nonlinear chain 
allows for the existence of rotobreathers (cf. Fig. lb). In 
the simplest case a rotobreather consists of one particle 
being in a rotating (whirling) state, while all others par- 
ticles with spatially decaying amplitudes from the center 
of DB: 

xoit + n) ^xo{t)+2Tr, 

To excite such a rotobreather the central particle needs to 
overcome the potential barrier generated by its two near- 
est neighbors. The rotobreather energy is thus bounded 
from below by Et, > 4. At sufficiently large energies 
Eb ^ 4 the central particle will perform nearly free ro- 
tations xo(t) ilbt + 5{t) with 5{t + Tb) = 5{t) being 
a small correction. The time-averaged off-diagonal hop- 
ping terms between site n = and n = ±1 are then given 
by 

(cos(fo(i) - xi{t))) « -{sm{Qbt){S{t) - xi{t))) (39) 

where xi{t) is also a small function. Thus the time- 
averaged scattering potential of a rotobreather presents 



a huge barrier that cuts the chain into nearly nonintcr- 
acting parts, contrary to the previous case of an acous- 
tic breather, where the acoustic breather potential cor- 
responds to a potential well. Therefore it becomes ex- 
tremely difficult for waves to penetrate across the roto- 
breather. The rotobreather solution with Qb = 4.5 is 
shown in Fig. 9. 
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FIG. 9. Displacements and velocities of an acoustic roto- 
breather at t = with fif, = 4.5. 

Inset: Relative strength for the first and second closed 
channels versus q. 

The corresponding Fourier components of the scatter- 
ing potential are shown in Fig. 10. 




FIG. 10. Fourier components Wnk for different k versus n 
for the rotobreather in Fig.9. 

In Fig. 11 we show the g-dependence of the transmis- 
sion for a rotobreather with Ob = 4.5 and compare it 
to the corresponding curve of an acoustic breather from 
the preceding section. We observe a dramatic decrease 
of the transmission at all values for the rotobreather, 
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in agreement with the above analysis. 

To decide whether the time-dependence of the roto- 
breather scattering potential is important or not, we first 
make a more precise computation for (39) (details of the 
calculations are presented in Appendix B) 

(cos(^o(i)-*i(*)))~-2^- (40) 

The time-averaged scattering potential for large fib cor- 
responds to two neighboring weak links of strength (40) 
inserted in a linear acoustic chain. Tlirec^ matrices Mn 
are enough to compute the transmission. Two of them 
involve matrix elements which are proportional to Cl^. 
The elements of the product will thus contain terms pro- 
portional to f2^, and according to (31) the result is 

iq ~ n-"". (41) 

This is precisely what wc also find from a numerical eval- 
uation of the transmission for the time-averaged scatter- 
ing potential in Fig. 12. However, although the transmis- 
sion coefficient tq for the full time-dependent scattering 
potential also drastically decreases with q and Qb, the 
dependence is weaker than for the case of time-average 
rotobreather potential. It scales as 

tq ~ n^^ (42) 

(see also Fig. 12). The reason is that besides a weak static 
link the rotobreather scattering potential has an ac term 
at frequency ilb of amplitude one (see Fig. 10). Thus an 
alternative route for the wave is to approach the roto- 
breather, to be excited into the first closed channel, to 
pass the breather and to relax back into the open channel. 
The corresponding scattering process of "virtual" absorp- 
tion and emission of phonons from the rotobreather can 
be also represented by three matrices as it happens for 
the dc analysis. However now instead of two weak links 
we have links of order one with a frequency change at 
site n = from Uq to iVq + Qb- This occurs in exactly one 
of the three matrices. Consequently the product matrix 
will contain elements proportional to and the trans- 
mission will scale as (42). 
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FIG. 11. The dependence of the transmission coefficient on 
q for Qt = 4.5. Data for acoustic rotobreather and aeoustie 
breather are shown correspondingly by solid and dashed lines. 




FIG. 12. The dependence of the transmission coefficient on 
the frequency of the breather Q), at fixed q = 0.21. Scattering 
by "real" rotobreather and a static part of DB are shown 
correspondingly by solid and dashed lines. 

We conclude this section by stressing that the wave 

scattering by an acoustic rotobreather is essentially re- 
lying on the time- dependence of the scattering potential. 
The rotobreather effectively cuts the chain in weakly in- 
teracting parts and thus hinders waves from propagation 
in a very strong way. 



V. OPTICAL BREATHERS 



In this Section we consider systems with a nonvan- 
ishing on-site potential V[xn] ^ (see Fig. Ic). The 
difference of such systems to acoustic models is the exis- 
tence of a gap in the spectrum of phonons |cjq=o| = I^"[0]. 
As a consequence the total mechanical momentum is not 
conserved, and the transmission coefficient now vanishes 
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not only at q = n but also at g = 0. An exception 
is the case when a localized Floquet eigenstate bifur- 
cates from the corresponding band edge for some special 
parameters [28 30]. Because of the presence of a gap 
in the plane wave spectrum there are now two different 
cases of interest - the breather frequency being located 
outside the spectrum > max|a;g| or inside the gap 
l^bl < min \LVq\. 

A. The case \Qb\ > max \uq\. 
Here we choose V{y) = ^y^ + ^y^ + jV^ and W{y) = 
In this case the spectrum of phonons is 

= l+4csin2(|) . (43) 

The breather profile for Clf, = 1.5 and c = 0.05 is shown 
in Fig.l3. 
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FIG. 13. The initial displacements of an optical breather 
with Qf, = 1.5 and c = 0.05 (velocities are zero). 
Inset: Sfc=i,2 versus q. 

The corresponding Fourier components of the scatter- 
ing potential are plotted in Fig. 14. 




FIG. 14. Fourier components Vnk for different k versus n 
for the breather in Fig. 13. 

We note that the DB scattering potential perturbs the 

diagonal terms and in this particular case {Qi, > ujq) 
presents a barrier for propagating phonons. For large 
breather frequencies the breather is strongly localized, 
i.e. essentially only one central oscillator is excited. The 
time-averaged scattering potential then becomes a single 
diagonal defect with a large strength /3o — f^l- It is 
straightforward to observe that the transmission coeffi- 
cient will thus scale as (see Eq.(25)) 

iq ~ n^^. (44) 

Due to the fact that the transmission vanishes exactly for 
both q = and g = tt we conclude that for large breather 
frequencies transmission is suppressed in general. 

However, in the case of an optical DB the time- 
dependent part of the DB scattering potential (more pre- 
cisely its second harmonic, see Fig. 14 ) is also large as 
f32 — for large values of the breather frequency. 

Thus, the time-dependent part of the DB potential can 
be rather important. Indeed, for the breather from Fig. 13 
the obtained transmission as a function of q shows that 
iq (see Fig. 15, dashed line) is at least one order of mag- 
nitude larger than tg (see Fig. 15, solid line). In addition 
tq shows a resonant minimum around q ~ 2.1 (sec in- 
set). Let us use the estimation (16) for k = 1,2. The 
parameters are wq.o = 3.5 , uq,! = 0.69 , t;o,2 = 3. First 
we find that si w 0.05 for all q implying that the k = 1 
closed channel does not participate in the transmission 
process. At the same time S2 > 10 for all q and thus, 
the k = 2 closed channel strongly participates in the 
transmission process. It is because the frequency of the 
propagating phonon 2Q{, — ~ 1.92 is close to a local- 
ized eigenmode of the time-averaged scattering potential 
Wl = 1-9. Notice here that the frequency of this local 
mode is above the phonon band. Thus we interprete the 
suppression of transmission as a strong coupling between 
the propagating wave and a particular localized mode of 
the time-averaged scattering potential, mediated by the 
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ac terms of the scattering potential (the k = 2 channel in 
this case). However, the simple estimation based on Eq. 
(16) does not show the resonant suppression of transmis- 
sion for a particular value of q and therefore, in order 
to quantitatively describe this effect we next apply the 
Green function method. 
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FIG. 15. The dependence of the transmission coefficient tq 
on the wave number q. The optical breather frequency is 
fib = 1.5 and the coupling c = 0.05. The results are shown 
for the time-averaged DB scattering potential (dotted line); 
the time averaged part and the second harmonic of the DB 
scattering potential (dashed line); the full DB scattering po- 
tential (solid line). 

Inset: same with tq on a logarithmic scale. Note the resonant 
suppression around q = 2.1. 

By making use of the estimation (16) we take into ac- 
count in (21) the terms with fc = and k = ±2, i. c. 
the scattering on the static potential and the phonon 
interaction with the second harmonic of the frequency 
2fli,. Moreover, we use the fact that the DB scattering 
potential is a strongly localized one, and obtain (the cor- 
responding diagram is shown in Fig. 3b) 

Ga;,(ni,n2) = G^^{ni,n2)+ 

+/32G^^(ni,0)GL™^(0,0)G.,(0,n2) , (45) 

where Go,, (ni , n2) is determined by Go, {rii, n^) that is the 
Green function of propagating phonons in the presence of 
time-average DB potential (the corresponding diagram is 
shown in Fig. 3a): 

G^(ni,n2) = G°(ni,n2)-^oG°(ni,0)Ga,(0,n2) . 



Moreover, the central part of Eq. 
G^Q^^^(0, 0), has a resonant form: 



1 



(46) 

(45), namely 
(47) 



where the local phonon mode frequency lo^ is mostly de- 
termined by a time-average scattering potential. How- 
ever, in the case as the phonon band is narrow (c ^ 1) 
the rcnormalization of lu^ due to the ac nonrcsonant pro- 
cesses has to be taken into account (the corresponding 
diagram is shown in Fig. 3c). The Eqs. (45) and (46) 
can be solved and we obtain 



Go,, (0,712) = 



G° (0,n2) 



l + (/3o-/3iGL^^U(0,0))GO^(0,0) 



(48) 



Thus, we arrive at the expression (25) for the transmis- 
sion coefficient tg but wit 
dependent parameter f3 



sion coefficient tq but with the renormalized wave number 



ujI - (206 - Wg)2 



(49) 



As the breather frequency is large (Oft » 1) both trans- 
mission coefficients for the static DB scattering potential 
and for the full time-periodic DB scattering potential, 
decrease with the breather frequency according to (44). 
However, the effective potential strength \/3\ is larger than 
Po and correspondingly the transmission tq is smaller 
compared to the transmission on the static DB scattering 
potential. Indeed, we observed this behaviour by direct 
numerical simulations (see Fig. 16). 

A most peculiar effect is that the presence of ac term 
in a scattering potential allows to tune the frequency of 
a local mode and to obtain the resonant suppression of 
transmission. Indeed by taking into account the type of 
diagrams shown in Fig. 3c we obtain for the renormalized 
local mode frequency ojl 



~2 2 , 



PI 



(40b - w,)^ - u)l 



(50) 



This formula is valid in the limit /?2 < 2'\/20j, [38]. For 
the particular case of the DB with frequency Of, = 1.5 
and p2 = i'02/2 = 1.5 we find that the resonant value of 
5 ~ 2 is remarkably close to the resonant suppression of 
tq around q = 2.1. Thus a strong dependence of lol on 
the amplitude of the ac part of DB scattering potential, 
and therefore, on the breather frequency, allows easily to 
change a position of the resonance. 

Note here, that taking into account the static part and 
the second harmonic of the DB potential allows to ob- 
tain a good agreement with the direct numerical simu- 
lations of a scattering by the "full" DB (compare solid 
and dashed lines in Fig. 15). It is also interesting to 
mention that the static; part and the first harmonic of 
the DB potential (fc = 1) lead to an increase of the 
transmission coefficient compared to the scattering by 
the time-averaged DB potential. It is just due to the in- 
terplay between the strengths of different channels of the 
scattering potential (see Eq. (49)). 
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Thus, we conclude that the presence of a closed channel 
allowing absorption and emission of phonons around the 
center of the breather leads to strong interference effects 
which are of destructive nature in the given example of 
optical breather. 




FIG. 16. The dependence of the transmission coefficient tq 
on the optical breather frequency Slf, for a particular value of 
wave number q = 0.5. The dashed line corresponds to the 
scattering by a time-averaged DB scattering potential. 



B. The case \i}b\ < min|a;q| 

Here we choose the on-site potential in the form V{y) = 
^y^ — ^y^ (note here that the results do not change if the 
cubic term has a positive sign) and the interaction term 
W{y) = The breather profile for ^j, = 0.85 and 

c = 0.15 is shown in Fig. 17. 
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FIG. 17. The initial displacements of an optical breather 
with f2b = 0.85 and c = 0.15 (velocities are zero). 
Inset: Sk=i,2 versus q. 



The corresponding Fourier components of the scatter- 
ing potential are plotted in Fig. 18. 




FIG. 18. Fourier components Vnk for different k versus n 
for the breather in Fig. 17. 



In this case the time-averaged DB scattering potential 
corresponds to a potential well in the diagonal terms, i. 
e. the parameter /3o in Eq. (46) has a negative sign. 
It leads again to the possibility to obtain a resonance 
as the frequency of an active closed channel matches 
the frequency of a localized phonon mode located in the 
phonon gap. Moreover, at variance to the optical DB 
with large frequency (preceding subsection), the num- 
ber of active closed channels may now increase substan- 
tially. This leads to a more complicated interference sce- 
nario between the open channel and several closed chan- 
nels. Indeed, the estimation of Sfe=i,2 shows that the first 
channel is strongly contributing, and the second closed 
channel can not be neglected either. At variance to the 
previous case for all propagating phonon frequencies the 
k = 1 closed channel has a much stronger contribution 
than the k = 2 one. However, we expect a resonant cou- 
pling between the propagating wave and the local mode 
through the k = 2 closed channel, as \u}q=o — '^^b\ = 0.7 is 
close to a local eigenmode of the time-averaged scatter- 
ing potential with frequency col = ±0.78. Moreover, the 
nonresonant processes in a strong k = 1 channel allow to 
renormalizc ujl similarly to the previous case of an optical 
breather with a large frequency. This resonant effect can 
be analyzed by making use of Eqs. (25), (45)-(49). In- 
deed, due to the resonance with a localized phonon mode 
the Green function G^^q^_,_^ (0, 0)^1 for a particular 
wave number qq. In this case the parameter (3 goes to 
infinity as the resonant condition, (2^^ — Wg^)^ — w£ = 0, 
is valid and correspondingly the transmission coefRcient 
vanishes. [39]. As the phonon frequency deviates from the 
resonant condition, the parameter /3 decreases and cor- 
respondingly the transmission coefficient reaches a max- 
imum (see Fig. 19). 
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In Fig. 19 wc show that the time- averaged DB scatter- 
ing potential provides with perfect transmission at a cer- 
tain wave number due to the presence of a quasi-bound 
state in the static scattering potential. At the same time 
the transmission for the full dynamical problem shows a 
maximum value of 0.1, and an additional minimum in 
t(q) with actually a full vanishing of transmission. These 
patterns are entirely absent in the scattering by the time- 
averaged potential. Thus we again conclude that the 
presence of active closed channels inside the breather core 
leads to strong interference effects. 




Q 

FIG. 19. The numerically calculated dependence of the 
transmission coofRcient tq (solid line) on q for the breather 
frequency Qb = 0.85 and c = 0.15. The dashed line shows 
the scattering by the time-averaged DB scattering potential. 
Inset: The same on a logarithmic scale. 

VI. DISCUSSION 

In this paper we considered the propagation of small- 
amplitude phonons through a nonlinear lattice in the 
presence of various discrete breathers. In particular, we 
obtained the transmission coefficient tq for the cases of 
acoustic and optical discrete breathers, and acoustic ro- 
tobreathers. We have shown that the presence of a DB 
leads to an effective scattering potential for phonons and 
this potential contains both static (time-averaged) parts 
and time-dependent periodic parts where the first and 
second harmonics are most important. Moreover, the 
strength and the width of the scattering potential are 
determined by the breather frequency f2{,. 

The static part of the effective DB scattering potential 
may fully describe the scattering outcome if the closed 
channels weakly contribute (acoustic breather with /3 = 
0). Equally the leading contribution to a finite trans- 
mission may come from a closed channel (acoustic roto- 
breather). A very interesting case is realized when a local 
mode of the time-averaged scattering potential resonates 



with a closed channel (acoustic breather with /? 7^ 0, op- 
tical breathers). This resonance leads to a full suppres- 
sion of transmission at some value of q. The nonresonant 
closed channels renormalize the local mode frequency and 
the corresponding position of the transmission zero in q. 
This effect has been discussed in several papers in relation 
to the Fano resonance [16,40]. A detailed explanation of 
the similarities and differences to the physics of a Fano 
resonances is beyond the scope of this work and will be 
published separately. It is worthwhile to note here that 
the location of a zero transmission value in q is not re- 
lated to the presence of so-called quasi-bound Floquet 
states, i.e. states which by parameter tuning are trans- 
formed from a localized state into one colliding (inter- 
acting) with the continous part of the Floquet spectrum. 
Instead a formally exact definition of a zero transmission 
is given in [28] through the asymptotic phase properties 
of Floquet eigenstates far from the breather center. 

The phonon scattering by inhomogeneities can be rel- 
evant with respect to the ongoing discussion about finite 
vs. infinite heat conductivity in acoustic chains with con- 
servation of total mechanical momentum Ref. [7-9,41]. 
The heat conductivity k mediated by noninteracting 
phonons is determined as 

K L j dqTg{L) , (51) 

where L is the size of the system. Here Tq{L) denotes the 
transmission of the given system. Thus, obviously in the 
absence of inhomogeneities the phonon heat conductiv- 
ity goes to infinity (k ~ L) as the size of the system L 
increases. However, the situation becomes more complex 
in the presence of inhomogeneities. In the case, when the 
inhomogeneities arc randomly distributed along the sys- 
tem, the number of inhomogeneities = nL, where n is 
the concentration of defects, and the total transmission 
Tq{V) ~ tq with tq being the transmission through a 
given defect. We obtain that the behaviour of heat con- 
ductivity is determined by the values of q where tq is close 
to one. Thus, the problem of an infinite heat conductiv- 
ity in the limit of a large size L naturally appears in the 
systems preserving the total mechanical momentum. 

In particular we can apply this general consideration 
to the phonon propagation in the presence of DBs. In 
the most interesting case of acoustic rotobreathers we 
obtain that the transmission coefficient tq ~ 1 — aq^ 
in the limit of small wave numbers. The coefficient a is 
determined by the frequency fif, and the specific choice 
of the acoustic rotobreather. Thus we argue that the 
mere presence of acoustic rotobreathers does not lead to 
a finite heat conductivity and k is still divergent in the 
limit of large L as 
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Here we assume that the heat is carried by phonons, 
and did not take into account phonon-phonon interac- 
tions. Another conclusion drawn from our results for 
acoustic breathers versus acoustic rotobreathers is that 
the g-domain around q = where the transmission is 
close to one is by orders of magnitude smaller for acous- 
tic rotobreathers as compared to acoustic breathers. This 
region is responsible for the quasiballistic transport of en- 
ergy by large wavelength phonons in the hydrodynamic 
regime and for the appearance of anomalous heat conduc- 
tivity. Numerical simulations for systems with acoustic 
rotobreathers have to be performed on length and time 
scales which are thus also orders of magnitude larger than 
the corresponding simulations for standard FPU chains. 

Our analysis can be also applied to the electromagnetic 
wave propagation through various Josephson transmis- 
sion lines. Although these systems are intrinsically dis- 
sipative, the dissipation may be rather small, and the 
transmission tq is still determined by the presence of in- 
homogeneities like rotobreathers [21,22,25,26] or dynamic 
edge states [42]. Moreover, in this particular case of 
Josephson lattices the properties of DBs can be easily 
tuned experimentally by changing the external dc bias. 
Notice here that local phonon modes located on a ro- 
tobreather play an important role in various switching 
processes [26] . A study of an electromagnetic wave prop- 
agation may allow a direct observation of this mode. 

Finally, a similar analysis can be carried out for a 
Schrodingcr equation in the presence of artificially ap- 
plied time-periodic perturbations [38]. This case is im- 
portant for the description of electron transport through 
a quantum wire in the presence of external electromag- 
netic fields. 
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APPENDIX A: TIME-AVERAGED SCATTERING 
POTENTIAL FOR ACOUSTIC BREATHER 

We take into account the oscillations of two lattice sites 
in the center of a DB only. These oscillations evolve 
exactly in antiphase with large amplitudes. The effective 
Hamiltonian is written as 

E=§ + § + \{xi+xt + {x,-X2)') (Al) 

and with xi = —X2, we arrive at the single oscillator 
problem with energy 



E = x' + -x" 



(A2) 



The frequency of oscillation (the breather frequency) is 
given by 



1 _ 2 /■ d 



M E ~ 1^* 



(A3) 



where Xm = \/^- After integration we obtain 



(A4) 



where B{x,y) is the S-function [43]. Next we compute 
the value oi X = ^3(a;^). We can express {x^) in terms 
of the energy 



{x^) = 



X dx 



(A5) 



After some algebra we find {x"^) 



9 B(|7|) ^"^^ 



x = ^|—nb■ 



(A6) 



APPENDIX B: TIME-AVERAGED SCATTERING 
POTENTIAL FOR ACOUSTIC ROTOBREATHER 

In order to calculate the time-averaged off-diagonal 
hopping terms we will take into account the central site 
which is in a rotational state, denoted by 4>, and two near- 
est neighbor oscillators (denoted by ai and a2). Thus we 
conserve the total mechanical momentum. The energy of 
such a system is 



£; = y -h 2(1 - cos((/) - a)). 

The two equations of motion are given by 

^ = — 2sin(^ — a), 
a = sm{<p — a). 



Introducing now variables u ■ 
rewrite the system (B2) as 

u = 0, 

w = —3 sin w. 
The relevant energy part is 



w 



I -|- 2a and w 



E = — cosw). 



(Bl) 



(B2) 



■ a we 



(B3) 



(B4) 



15 



Now we compute the average coupling between rotational 
and oscillatory states for large frequencies 



e = (cosw) = 




cos w dw 



■COS It;) 



(B5) 



Because £^ « we find 



271- 




/3(1 — cos to) , , , 

(1 + = -) cos w dw (B6) 
2E ^ 



which leads to 



and finally to 



4y2£;3/2 



(B7) 



(B8) 
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